Kinetic surface roughening for the Mulhns-Herring equation 
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Using the linearity property of the Mulhns-Herring equation when the velocity is zero with a 
Gaussian noise, we obtain an analytic form for the global mean-square surface width and height- 
height correlation function. This can be used to read the critical exponents in any dimension. In 
particular for d = 1 we show that although the surface is super rough the system exhibits Family- 
Vicsek scaling behavior. 
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I. INTRODUCTION 



_ In nature one can find a huge number of systems that develop a rough interface in the process of growing. Many 

I of them have been understood by the use of some tools from fractal geometry [J , such as scaling analysis as well as 
^ through modeling with stochastic partial differential equations (SPDEs) and discrete models. While these concepts 
I t ido not constitute an exclusive theoretical framework for surfaces in the physical world, they can be applied to gain a 
Ci 'deeper understanding of many processes, for example, in biology. Due to the many possible important applications 
^ ,in medicine, tumor growth constitutes one of the most interesting subjects of study to which scaling analysis can 
be applied. Actually very important research on tumor growth has been recently carried out. Strong experimental 
evidence has indicated that tumor growth belongs to the molecular beam epitaxy (MBE) universality class 0, . 

The MBE dynamics is characterized by a number of features that include a linear growth rate, the constraint of 
growth activity to the outer border of tumor, and surface diffusion at the growing edge; all of them have been observed 
experimentally. Stochastic models for tumoral growth have recently been considered in [1, Q . 

The continuum equation which describes the MBE universality class, known as the Mulhns-Herring equation 
> has the following SPDE: 
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^ = -nV'h + Fo+vi^,t), (1) 

where h is the interface height, and k is the surface diffusion coefficient. Fq has dimension of velocity and may be 
interpreted as the product of the mean cell radius and the cell division rate. 7/(x, t) is a Gaussian noise with zero 
mean and 

(77(x, t)r,i^',t')) = 2DS^{^ - ^')S{t - t'). (2) 
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'The critical exponents can be extracted from the equation ([T]) simply by power counting. For d = 1, if we ignore for 
the moment the velocity Fq, and perform the transformations x — > bx, t — > h^t, and h — > h°'h, we will get the scale 
'invariant system only ii z — A and a = 3/2, indicating that the system is super rough. Actually in this particular 
case the system is characterized by a set of critical exponents given by a = 3/2, aioc = 1, z = 4, /3 = 3/8, and 
/3* = 1/8 Q • Interestingly enough it was found that this system with the above critical exponents is compatible with 
the experimental data obtained from tumor growth measurement 0, Q . 

One nice feature of the Mulhns-Herring equation is its simplicity; it is linear and can be solved exactly by means 
of Fourier transformation. 

We will be interested in the dynamics of a d dimensional surface, which will be described by the stochastic variable 
h{x,t) giving the height of the surface at time t above substrate position x. Starting from a initially flat substrate of 
linear size L at time t = 0, the global mean square width of the surface. 



with N being the total number of substrate sites, satisfies [8| 

w^iL,t)r. L'-f(±]. (4) 
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For t <^ the L dependence drops out, and we get w^{L, t) ~ t^^ . On the other hand for t ^ the scahng function 
becomes constant, and the Eq. ^ yields w'^ ^ L^" Q. 

Another interesting quantity we wiU study is the height-height correlation function 

C{r,t) = (^\h{r,t) - h{0,t)\^) , (5) 

which scales in the same way as the square of the local width, in the form of 

2a J ' 



Cir,t)^r^- U^-j (6) 

where r = |r| and the scaling function fd{u) behaves as 

{const if M 3> 1 

(7) 
u'^'^ if u <C 1 

Studying the asymptotic behavior of the height-height correlation function provides an alternative method to deter- 
mine the critical exponents. 

The aim of this article is to calculate the mean square surface width and the scaling form for MuUins-Herring 
equation. We will also find an exact expression for the height-height correlation function C(r,t) with an emphasis 
on the asymptotic behavior of the correlation function. Using this expression we find the scaling function fd{^) for 
different dimensions d. In particular for d = I where the model could be used to describe the tumor growth we show 
that the surface is super rough and exhibits the Family- Vicsek behavior. 

II. THE LINEAR LANGEVIN EQUATION 

Let us start from the linear Langevin equation which can be read from ([1]) by setting Fq — 

rlh 

^ = -nV'h + fj{^,t), (8) 



Although the exponents 



2 = 4 (9) 
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which characterizes the dynamics scaling of this equation are well-known, it is interesting to determine the value of 
various quantities in the scaling regime. These values are related to each other through the parameters n and D of 
the continuum model. 

Using the Fourier transformation one can find a solution for Eq.® as follows 

h{k,u;)^ V 7?(k,c.), (10) 

where fc = |k| and 

/i(k,w) = j d'^x dt e^^e^-^-^'^/iCx,*), 

r]{k,uj) = J d'^x dt e-*(''-^-'^*)77(x,t). (11) 

We consider a growth process starting from a flat substrate in which one should set ?7(x, i) = for t < 0. The noise 
correlator in Fourier space then takes the form 

{r]{-k,ij)Tj{k',uj')) = 2D{2TTf5\k + k') dTe'^^+^'^^ (12) 



By making use of the Eqs. (fTT|) and (fT2|) and after performing the w integrals, we get the following expression for the 
amplitude /i(k, t) 

(/i(k,t)/i(k',t')) = -?7(l-e-2-fe't)(27r)'',5''(k + k') for t>Q (13) 
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A. Mean square width 



In this subsection we obtain a closed form for the mean square width which describes the interface fluctuations and 
study its asymptotic behaviors. This could be used to read the critical exponents. 

To proceed we note that the mean-square width ([3]) can be recast to the following form 



w^{L,t) 



{h{k,t)hik',t)). 



(14) 



Performing the momentum integration over a spherical shell from ^ to ^ and carrying out the angular integrations, 
we arrive at 



dk D 



(15) 



where K^'^ = 2"^-^ ^r^rd). In comparison with the lattice model, a should be identified with the lattice spacing. 
Using integrating by parts, from Eq. (|15p one finds 



D 



Kt , 



where 



and 
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d - A \ K I 



Fdix) 



4-d 



(4-d)(2^) 



^ , , r(l ) 1 ■ (8-5d) d (4-d) 



-1 ^-y 



dy 



1- e(-32-M + (327r4x)i-3 
For large a;, Fd{x) approaches a constant, while at small x, the limiting form of Eq. (|18[) is given by, 

0{x) for d > 



r(|) 4-d 



(16) 



(17) 



(18) 



(19) 



leading to the exponents as in Eq. ([9]). 

We note, however, that for d = 4 special care is needed. Actually since in this case a and (3 are zero one needs to 
consider this case separately and indeed redo the computations from first step with d = 4. It is then straightforward 
to go through the computations to find 
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In 



L 



F 
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where 



In F{x) = - \si2 + -Ei (-32 tt^ x) 



(20) 



(21) 



Here Ei{—y) = — u ^ e " du is the exponential integral. By making use of the asymptotic behavior of Ei{~y) 
for small y, one gets 



w^{L,t) 



D 



7 + In 



27r*'Kt 



for a < [Kt) < L 



for [nt) " > L 



(22) 



In this case at early times the width scales logarithmically with time and the saturation width depends on the 
logarithm of the system size. 
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B. Height- height correlation function 



The aim of this subsection is to reproduce the critical exponents using the height-height correlation function. To 
do this we will first obtain a closed form for the height-height correlation function, then we will read the form of the 
scaling function. By making use of the asymptotic behavior of the scaling function we can read the critical exponents. 

To proceed we note that by going to the momentum space and using Eq. ([Tn|) the height-height correlation ([5]) 
reads 



C(r,t) 



d'^k 2D 

Carrying out the angular integration, one finds 



1 



1 — COS (k.r) 



C(r, t) - Kd 



■K j a 



2tt/L 



dk k 



-2Kfe*t 



Jd/2-1 (kr) 



(23) 



(24) 



where r — |r|, and Jn{x) is the nth-order Bessel function. In general we need to evaluate this integral and study its 
asymptotic behaviors as a function of r. To do this we note that the integral can be broken into two parts. The first 
part is independent of r, and the r dependence comes from the second part. More explicitly, assuming r <^ L and 
(k^)^/'' <^ L, one has 



2n 

C{r,t)=Kd — 

K 



7r(Kt)i/*/a 



dz 



„4-d 



7rr/ a 



dyy (2/,)i-^/2 [} 



(25) 

Here we have used the change of variables z — (Kty^^k and y = rk. The first integral can be carried out leading to 



2D 



{Kt)- 



[d - 4)7r4-'^ 2i+'^/4 



r(d/4 - 1) - r(d/4 - 1, 2T:^Kt/a^ 



K 



r nt 



(26) 



for d ^ A. Here Qd is defined by the second integral in (|25p . For a <^ r this function is independent of a/r and the 
height-height correlation function can be recast to the following form 



2D 



C(r,i) = — r^--*/, - +0(1), 



Hit 



(27) 



where fdint/r'^) is the scaling function. 

It is then easy to evaluate the scaling behavior of the height-height correlation function which for d 7^ 2 is given by 



rid — 



C(r,t) 



id- 5 



(d-4)7r<i-<* 



T{d/2)T{d/2 - 2) r 



4-d 



2D 



r(d/4-l) 

2i+ti/<i 



{Kty 



for a < r < {kI) ^ 
for a < [nty < r 



(28) 



For d = 2 the calculation yields 



f4,(0(l) + log(^)) 



C(r,t) 



^ i^tf" ( 0[l) - O 



t\ll/6 



for a 4; r <C {nt) 4 
for a <^ {Kt)i ^ r 



(29) 



We note that there are two extra terms appearing in the asymptotic behaviors (logarithmic term for -pi :i> 1 and 
exponential growth for ^ <C 1) which make it impossible to find the scaling form for this case. 

We note also that in the case of d = 5 the roughness exponent is negative indicating that the interface is flat. 
Taking into account that for the case of d = 3 the exponent is positive one may conclude that there would be a phase 
transition for d = 4 case. Indeed this is exactly what we have observed from mean square width's computations. In 
fact for d = 4 the Eq. ([24]) reads 



C(r,i) 



1 D 

4 7r2 K 



Tr/a 



dk k 



-2.k- t\ f ^Mkr) 
\ kr 



(30) 



This integral can be evaluated leading to the following limiting behavior 

f 4^ f log ^ + ^(1) for a « r « iKt)i 

C{r,t) = \ (31) 
I 16^ f log(fj) + Oil) for a « {Kt)i « r 

As we see the correlation function decays logarithmically and therefore we get the logarithmic behavior for the scaling 
function, as expected. 

III. CONCLUSION 

In this paper we have obtained explicit expressions for the mean square width of the d dimensional surface for the 
Mullins-Herring equation which in asymptotic behavior has the standard Family- Vicsek scaling form. We have also 
considered the height-height correlation function and studied its asymptotic behavior to read the critical exponents. 

Although we have worked out the height-height correlation function for any dimension, the most interesting case is 
d = 1 where from ([28|) we find a — 3/2, z — A and /3 = 3/8. Actually it was shown in 0, Q that these exponents are 
compatible with the experimental data obtained from tumor growth measurement. In other words one would expect 
that the one dimensional MuUin-Herreing equation can describe the tumor growth. Therefore it can be utilized to 
study the roughness of the tumor surface using, for example, the height-height correlation function. In particular we 
have seen that since the global roughness exponent is bigger than one the surface is super rough and the height-height 
correlation saturates. This shows that the surface of tumor is super rough and exhibits the Family- Vicsek behavior. 
This has to be compared with the numerical results presented in 7] it was shown that for large L the system does 
not display any anomaly indicating that scaling function has the Family-vicsek shape. 

Finally, for completeness, let us also consider the auto height correlation in the limit t, t' —>■ oo, 

{\h[v,t) -h{v,t'f'^ = Kd rffcfc'^-^ ^ X (l-e-"'^' (32) 

Performing the steps similar to the previous case we get 



16 K 



^log(^^F^) + for d = 4 



(33) 



where A is given by Eq. ([TT]). 
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